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We consider a general approach to the theory of continuous media 
starting from Lagrangian formalism. This formalism which uses the 
trajectories if constituents of media is very convenient for taking into 
account different types of interaction between particles typical for 
different media. Building the Hamiltonian formalism we discuss some 
issues which is not very well known, such as relation of famous Thompson 
theorem with the symmetry with respect to volume preserving diffeomorphisms. 
We also discuss the relation between Euler and Lagrange description 
and present similar to Euler C 2 formulation of continuous mechanics. 
In these general frameworks we consider as examples the theory of 
plasma and gravitating gas. 

1 Introduction 

At the present time almost all fundamental physical phenomena could be 
formulated in the frameworks of either classical or quantum mechanics. That 
means that these phenomena admits the Hamiltonian description, which due 
to its long history developed many powerful methods of analysis of the general 
properties of evolution of the systems and the tools for the solutions of partial 
problems. 

In this respect the fluid mechanics stands aside (in spite of its name) from 
orthodox mechanics. The reasons for that is not only the infinite number of 
degrees of freedom of fluid which could be treated e.g. within classical field 
theory or statistical mechanics. The main difference between conventional 
field theory and continuous mechanics is that in the first case we can speak 
about the dynamics of the field at one point in space (which of course 
interacts with the field at the neighboring points), while in the case of the 
fluid, describing the interaction of the neighboring particles which constitute 
the fluid we are loosing its position in the space due to the motion of fluid. 
In the same time the objective of usual problems of hydrodynamics is to 
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define the velocity, density and a thermodynamical variable (pressure or 
entropy) as the functions of the coordinates x and time t for the appropriate 
boundary conditions and/or initial data [Tj. The similar problems also appear 
for magnetohydrodynamics dealing with sufficiently dense plasma [2j. For the 
developing of the Hamiltonian formalism we need to start with more detailed 
description based initially on the trajectories of constituents of media. This 
description is especially important for plasma, because the fundamental electromagnetic 
interaction could be formulated only in terms of the trajectories of the 
charges. Needless to say that some aspects of this approach was extensively 
studied in the series of papers by J.Marsden, A.Weinstein, P.Kupershmidt, 
D.Nolm , T.Ratiu and C.Levermore [5] especially in the context of stability 
problem. 

Of course not all properties of fluid could be formulated in the frameworks 
of the Hamiltonian approach. For example, we leave open the question of the 
energy dissipation, viscosity et cetera. 



2 The Lagrangian equations of motion 

In fluid mechanics there are two different pictures of description. The first, 
usually refereed as Eulerian, uses as the coordinates the space dependent 
fields of velocity, density and some thermodynamic variable. The second, 
Lagrangian description, uses the coordinates of the particles x(£i,t) labeled 
by the set of the parameters £ i5 which could be considered as the initial 
positions £ = x(t = 0) and time t. These initial positions £ as well, as 
the coordinates x(£i,t) belong to some domain D C i? 3 . In sequel we shall 
consider only conservative systems, where the paths of different particles do 
not cross, therefore it is clear that the functions 2?(&, t) define a diffeomorphism 
of D C R 3 and the inverse functions £(xi,t) should also exist. 



The density of the particles in space at time t is 

p{x, t)= [ d 3 ^p {^)5{x - *)), (2) 
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where Po(0 is the initial density at time t = 0. The velocity field v as a 
function of coordinates x and t is: 

v(xi,t) = x(£(xi,t),t), (3) 

where £(x,t) is the inverse function (pQ). The velocity also could be written 
in the following form: 

or 



t)v(xi, t) = J d £po(&)£(6, t)8{x - x(ii, t)). (5) 
Let us calculate the time derivative of the density using its definition ([2j) : 

/d 

= J d 3 Z Po &) ^5{x-x{i h t)) 

= -£/ d^p ^mu t)5{x - a?(&, t)) 

d 

= -—p{xi,t)v(xi,t) (6) 

In such a way we verify the continuity equation of fluid dynamics: 

p(xi, t) + d(p( Xi , t)v( Xi , t)) = 0. (7) 

Using the coordinates x(£i, t) as a configurational variables we can consider 
the simplest motion of the fluid described by the Lagrangian 

L = /^(«^P. (8) 
The equations of motion which follow from j8]) apparently are 

mpa(&Z(Zi,t)=0 (9) 

Now let us find what does this equation mean for the density and velocity 
of the fluid. For that we shall differentiate both sides of (jSJ) with respect to 
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time 



d 

dt' 



p(xi,t)v(xi,t) = j d 3 £p (£i)x(£i,t)5(x - x(£i,t)) 

d 

+ J d 3 tpo(ti)x(£i,t)—5(x - £(&,<)) 



(10) 



The first term in the r.h.s. of (fTOl) vanishes due to the equations of motionQ 
and transforming the second in the same way, as we did in ([6]) we arrive at 



d d 

—p(xi,t)v(xi,t) + —(p(xi,t)v(xi,t)v k (x h t)) = 

Let us rewrite (TTT1) in the following form: 





0. 



(11) 



v(Xi,t) p(Xi,t) + -—(p(Xi,t)v k (Xi,t) 

1 ox k v 

+p{Xi, t) \v(Xi, t) + v k (x h t)—V{Xi, t) 



(12) 



The first term in (1121) vanishes due to the continuity equation, while the 
second gives Euler's equation in the case of the free flow: 



d 

v(xi, t) + V k (Xi, t)- — V{Xi, t) = 
dx k 



(13) 



Before moving further we need to make one comment concerning the 
Lagrangian (JSj) . The reader may get an impression that this Lagrangian 
depends not only on the dynamical variables at the time t, but also at the 
initial time t = through the presence of po(6)- I n order to eliminate the 
doubt let us take the unity 
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d s x5(x — x(£i, t)) 



(14) 



and insert it into the integrand of 



d s £d 3 xp (£i 



mx (d,t) 



8(x- Mis- 



performing the integration in (1151 ) over £ we obtain: 

d*xp(x- t ) ™ {Xi ' t)2 



(15) 



(16) 



This representation of the Lagrangian (JSj) shows explicitly the dependence of 
L only of the dynamical variables at the time t (besides we are not going to 
consider here the Euler variables as an independent). 

In order to have more realistic model of the fluid we need to introduce 
into the Lagrangian ([HI) the "potential energy"term which will give rise to 
the internal pressure field in Euler's equation. 

As we have mentioned above, the functions t) define a diffeomorphism 
in R 3 therefore the matrix 

4(6.*) = ^#^ (17) 

is non-degenerate for any and t. The integral (pQ) may be expressed via the 
Jacobean - the determinant of A 3 k (£i,t): 



p(x t) Pote) 



detAfat) 



(18) 



For some problems one can choose the initial density Po(£i) to be uniform 
in D C R 3 1 and effectively normalize the density field p(£i,t) by putting 
Po(Ci) = 1 ( one particle in the elementary volume). However, we prefer to 
keep Po(£«) arbitrary, because in the case e.g. the presence of soliton in the 
fluid the density could not be uniform at any time. 

Taking the "potential energy"to be the functional of detA we can write 
the lagrangian in the following form: 



mf - f{detA{^t)) 



(19) 



In order the Lagrangian be the functional of the dynamical variables at the 
time t, the function f(detA(^ t)), should has the following form: 

/™> = n^)^ (20, 



Now the equations of motion become 

mpoi&xjfat) - J-(po(&) (A' 1 ). k {i u t)f\detA)detA \ = (21) 
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Substituting x'j(^,t) from ( 11 39[) to the equation (jTOl) and acting as we did in 
the derivation of the equation (fTBT ). we obtain 



( d d \ d 

mp(x h t) ( — + v k (xi,t)— j Vj(xi,t)-— [ pQ(£i)f(detA(£i,t)) 




It is now obvious that if we identify the — ^p (^i) /'(detA^i, £)) _ , 

with pressure p(xi,t), the equation (|140l) takes the form of usual Euler 
equation without viscosity: 

p(xi,t) + v k (xi,t)-^j Vj(xi,t) = —^-p(xi,t). (23) 

For the function f(detA(£i,t)) parameterized as in ( |20l ). the pressure has the 
following form: 

The pressure p(xj, t) which appeared here is the result of interaction between 
particles, which constitute the fluid. We can also add to the Lagrangian (fT9l) 
the term, which describes the interaction with an external field: 

L ext = ~ J (PtPai&Uixifo, t)), (25) 
and the Euler equations take the form: 

^(s^-O^-^-^^- 1261 

As we shall see later the description in terms of Lagrange variables is 
very convenient for introducing interaction between particles because in this 
picture the usual coordinates are the dynamical variables. For example, 
considering two components fluid - ions and electron, interacting with electromagnetic 
field we obtain the theory of plasma. Introduction the interaction of particles, 
which constitute the fluid with color gluon field gives us the theory of quark- 
gluon plasma. Very interesting system is the gravitating gas, where particles 
interact with each other through gravitation field. These different types of 
interaction changes the physical properties if the media, but all of them could 
be described within the common approach. 
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3 Hamiltonian formalism 



The description of fluid in terms of Lagrangian variables brings no difficulties 
in construction of canonical formalism. Indeed, for the Lagrangian ([19]) . 
the canonical coordinates will be the functions x(£i,t) and its conjugated 
momenta are defined as the derivatives of the Lagrangian with respect to the 
velocities x(£i, t): 

5L 

P&' £ ) = * = mpoi^M^t), (27) 

The Hamiltonian is given by the Legendre transformation of the Lagrangian: 



H 



< 28 » 

The canonical Poisson brackets is defined by 

tete),^(4')} = M 3 (6-0 (29) 

Apparently the Poisson brackets (|29l) and the Hamiltonian (l28l) define the 
equations of motion for the canonical variables, which are equivalent to the 
Lagrange equations. The phase space Y of the fluid is formed by t),p(£i, t). 
What is missing at the moment is the x-space interpretation of the variables 
x(£i, ) and p(£i) and the goal of this section is to find a canonical transformation 
of x(£i, ) and p(£i) to equivalent set of coordinates in T, which are x-dependent 
functions. 

Let us introduce the new objects using the same averaging, as we used in 
the previous section 

l(x) = [ d^p^x-x^)) = p{x)p[^{x)). (30) 



The Poisson brackets of l(x), induced by (1291) has the following form: 

{lj(xi),lk(yi)}= [ l k(xi)^ + lj(yi)-^\5(x - y). (31) 

These Poisson brackets was introduced geometrically as "the hydrodynamic- 
type brackets"long ago in the papers [7] without physical explanation. The 
present discussion reveals the origin of these brackets. The commutation 



7 



relation (I3TT) coincides with algebra of 3-dimensional diffeomorphisms. In 
other words, l(x) are the generators of the finite diffeomorphism Xj — > <j>j{xi) 
of any x-dependent dynamical variable in Y. It should be mentioned that the 
group of diffeomorphisms, generated by lj(xi) is not a gauge symmetry in case 
of fluid mechanics, as it is in the case of e.g. relativistic string or membrane. 
In the same time in fluid mechanics there is an infinite dimensional symmetry 
(but not a gauge symmetry) with respect to special (i.e. volume preserving) 
diffeomorphisms SDiff, which will be considered in later. 

Let us leave for a moment l(x) and consider another x-dependent functions 
£,i(x) which are the inverse to Xj(£i) functions (CTJ) . Differentiating the first 
equation (jTJ) with respect to x we obtain: 

A j k (Ux k ))^^ = 5j, (32) 
in other words, the matrix 

°fte) = ^r 1 (33) 

is inverse to A 3 k (£i(xk))- Therefore, from (11361) it follows that 

p(xi) = p {£i{xi))detai(xi) (34) 

To simplify the calculation of the Poisson brackets we can express £, (£j) in 
the following form: 

t / n _ J d 3 tpo{£i)tj5{x - x{£j,t)) 
From (j35|) we easily obtain 

{£i{xi)MVi)} = Q (36) 

The calculation of the Poisson brackets between £j(xi) and lj(xi) is more 
involved, but the result is simple: 

{ii(^),&(w)> = -^^ 5 ^-y) (37) 

In such a way we have constructed the set of x-dependent coordinates (lj(xi), Cj( x i)) 
in the phase space T, but the transformation 

M&,xA&)^(hfa)M*i)) (38) 
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is not canonical. The set of canonical ^-dependent coordinates in T could be 
obtained in the following way. Let us multiply both sides of ( 1371 ) by matrix 

Aim^mA^^kiVi)} = S k J(x-y), (39) 



where we have used (1321) . Further, due to the relation (1361) . we can put 
Mn^ii. x k)) inside the brackets and obtain: 

{7r m (x i ),Uyi)} = ^8(x-^, (40) 

where 

^m(Xi) = -A-^^Xk^ljiXi) (41) 

By tedious, but direct calculation we also obtain 

{n m (xi),7r k (yi)} = 0, (42) 

so the set (vr m (xj), ^k{Vi)) i s formed by the canonical variables. In terms of 
these canonical variables the generators of the group of diffiomorphisms lj(xi) 
has the following form: 

lj(Xi) = -^^-7r k ( Xi ). (43) 

For the Lagrangian considered, from ( 1271 ) and ( 1301) follows that lj(xi) is given 
I'.v 

lj(xi) = mp(xi)vj(xi) (44) 

and the representation (1431) is very similar to Clebsh parametrization [8] of 
the velocity. The distinction of (T431 from the original Clebsh parametrization 
is the appearance of three "potentials instead of two, for the 3-dimensional 
fluid. The reason of this difference will be discussed in the end of this section. 

The canonical Hamiltonian ( 1281 ) should now be expressed in the terms of 
new variables (7r m (xi), Ck(lJi))- Indeed, let us again insert the unity 



1 = / d s x5(x -£(&)) (45) 
into the integrand (|28|) and change the order of integration: 

H = Jd 3 xJ d 3 £ ^ —±^ f^ h t) + f(detAfo))\ 6(x-x(^)), (46) 
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Performing the integration over £ with the help of (j2j) and (IBTil ) we obtain: 

WM^^+^rfb 5 )' <47) 

Making use of (1331 ) we can express H via canonical variables {ir m (xi), £fe(?/i)): 

w - / A ( aJt^^y ^'-^ + ' <48) 

where the function V(p(xi)) for the "potential "part of the energy was introduced 
in ( |20l l. For usual fluid or gas this term represents the internal energy of the 
fluid and it should vanish for uniform density distribution p(xi) = p as . A 
phenomenological expression for V{p{xj)) could be written as follows |9J: 



V{p{x i ))^^{5p{x i )f + X{Vp{x i )f + ... , (49) 
where Sp(xi) is the deviation of the density from its homogeneous distribution: 

6p(Xi) = p(Xi) - Pas- (50) 



The first term in (1491 is responsible for the sound wave in the fluid (k is the 
velocity of sound), the second term in (1391 describes the dispersion of the 
sound waves. 

In order to reveal the relation of the Hamiltonian flow, generated by (1351) 
with geodesic flow [4j let us introduce the metric tensor gjk(xi): 

g jk (xi) = A^(x k ))A^(x k )). (51) 
With this notation (1351) takes the form: 



H= fd 3 x[ -^ V /^)/ m (x i ) 7 r fe (x i )7r m (^) + V{-^=) . (52) 

The metric tensor with upper indices g^ k (xi) denotes, as usually the inverse 
matrix and ^ 

g(xi) = detg jk (xi) = (53) 
fr[Xi) 



The representation (1521) permit us to consider the hydrodynamics as the 
geodesic flow on the dynamical manifold with metric gjk(%i) and many general 
properties of the hydrodynamics could be derived from this fact (see e.g. |4J). 
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4 Infinite-dimensional symmetry and integrals 
of motion. 

As we have mentioned above, the Lagrangian ( fT9l) possesses the invariance 
with respect to the "volume preserving"group of diffeomorphisms SDif f[D], 
where D C R 3 . Indeed, let us write the Lagrangian ([191 ) in terms of the 
Lagrangian density (to simplify the equations, in this section we shall take 
the initial density po(&) — 1 ) 

L = J d 3 ££(&) (54) 
and consider the transformations from SDiff[D] of the coordinates e D 

(55) 



det ^l^lL = i. ( 56 ) 
Apparently, due to ( 1561 ) we obtain: 

L = J d 3 ^ t )C (0(6)) = / d 3 £C (0(6)) (57) 

and according to Noether's theorem this invariance results in the existence of 
an infinite set of integrals of motion. To obtain these integrals we first need 
to find the parametrization of the transformations (|55l) . (|56l) in the vicinity 
of identity transformation: 

= + ^(6)- (58) 

From (l56l) follows the equation for (Xj(£i): 



0. (59) 



Further we must explicitly take into account that the volume preserving 
diffeomorphism (!55l) leaves the boundary of D invariant. We shall limit 
ourself with the case when D is formed by extraction of the domain with 
the differentiable boundary given by 

<?(&) = (60) 
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from R 3 . Physically that means that we put in the fluid the fixed body, 
the shape of which is given by (1601) . The condition that the infinitesimal 
diffeomorphism (1551) preserves D in this case is 

0& + a i(6))|, (6HJ = °. ( 61 ) 

or 

«,te)V^(6)| fl(Ci)=0 = 0. (62) 

Geometrically equation (l62l) means that the vector is tangent to the 

surface, defined by f l60l ). because the vector Vj#(£i) is proportional to 

the normal rij(£i) of the surface ( l60l) at the point £j. 

From Noether's theorem we obtain that the invariance of the Lagrangian 
with respect to the transformation ( 1551 ) gives the following conservation law: 



| / d 8 fen(&)^^a l (&) = 0, (63) 

where a{(&) satisfies to the conditions (1591) and ( |62l) . The existence of these 
conditions forbids to take the variation of the l.h.s. of (1631 ) over cti(£i) and 
obtain the local form of integrals of motion. For that we need to extract from 
(159]) and (l62l) the integral properties of (fy(£j). 

Consider an arbitrary, single- valued, differentiable in D function 
Then the following equations are valid: 



JD Otj 



The first equality is valid due to condition (l59l) . Using Stokes theorem we 
can transform the integral of total derivative in the last equality ( |64l) : 

J D d^{m^m) = f^ds^m^M) = o, (65) 

where 3D denotes the boundary of D. The last integral in (1651 ) vanishes due 
to condition (1621) because the differential dSj is proportional to the normal 
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vector of the surface dD, defined by ( l60l) . From (1641) and ( 1651) we conclude 
that 



^(6)^ = (66) 

for any smooth, differentiable Taking this property of «? (£;) into account 

we obtain from the conservation laws (1631 ) that the quantities 



Jk&) = Pm&)^^ (67) 

are conserved modulo some term which is the gradient of a scalar. In particular, 
that means that 

Rj&) = ^i-^-m) = e ife , J- (^(6)%^) ■ (68) 

is the integrals of motion. Note, that as the group of invariance is infinite- 
dimensional, we obtain an infinite number of integrals of motion. With respect 
to Poisson brackets f [29l the -Rj(£j)'s form an algebra. This algebra could be 
written in a compact form for integrated objects 

R[<f>] = J d 3 £&(&)^(6), (69) 

where (/)](&) are smooth, rapidly decreasing functions. The algebra of R[4>] 
induced by Poisson brackets (|29l) has the form: 

{R[(f)), R[ip)} = R[curl<f> x cw#] (70) 

The construction of the x-dependent object, corresponding to Rj(Ci) is n °t 
an easy task, because our "averaging "with 5(x — x(£i)) will introduce time 
dependence and instead of conserved object we shall obtain a density, whose 
time derivative gives a divergence of a "current". Therefore we need to 
introduce another kind of averaging without explicit refereing to the x(^) 
coordinates. For that recall that under diffeomorphism a closed loop transforms 
into closed loop. Then let us consider such a loop A and a surface a whose 
boundary is A. The integral 



V=[dS j R j &) (71) 



13 



where the vector dSj is as usually the area element times the vector, perpendicular 
to the surface, is conserved, because of the conservation of Rj(£i). Further, 
from Stokes theorem we have: 

r {^ jPm ^) d -^-. (72) 

Changing the variables in (1721 ) we obtain: 

T — / dxj J 'f ~r = / (ix.f j(xi), 

where A is the image of the loop A under diffeomorphism ^ — > Xj(£i). The 
object (1731) is very well known in hydrodynamics as the "circulation "and 
its conservation is known as W.Thompson theorem [10]. The relation of the 
circulation conservation with the invariance under special diffeomorphisms 
was first explicitly established in [H], though it also could be extracted from 
general discussion in Appendix 2 of [4]. In order to avoid confusion we should 
make a remark concerning the meaning of the contour integral ([TBI) . The 
matter is that the diffeomorphism ^ — > Xj(^i) transforms a fixed loop A into 
time dependent loop A and while e.g. calculating the Poisson brackets of T 
with Hamiltonian we need to differentiate not only the integrand, but also 
the contour. The explicit form of this integral is 

T = f A ds dXm ( | ( S) ' t] v m (x k (Us),t),t), (74) 

where £»(s) is the contour A. In this form it clearly seen that the circulation 
could not be expressed in terms only Euler variables. In next section we shall 
come back to this issue and discuss the relation of the symmetry with respect 
to SDiff and Euler description. 

Conservation of circulation is not the only consequence of ( 1631) . Consider 
for example the case of 2-dimensional space. Here, instead of the conserved 
vector Rj(£i) we shall have the conserved scalar 

= e« J-Ji(6) = e« £- (pU&^P) , (75) 
This scalar defines the following integrals of motion: 



(73) 
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Changing variables in ([76]) £j — >■ ^(xj) which is possible, because £,( IS el 
diffeomorphism of D, we obtain: 



1 n 



D d 2 ^o^Ai(e(x))Ar(e(x))^^j , (77) 



where the matrix was defined in (fT71) . We can present J n as 



/. = / d ^)-(^.A^)". (78) 

using the following property of 2-dimensional matrix A 3 k (£): 

e kl A{(£(x))A?(ax)) = e jm detA{^{x)) = e jm -^, (79) 

Using the relation between l(xi) and velocity we can rewrite the integrals of 
motion in case of 2-dimensional fluid in the following form: 



d^pix^idM^) ~ d 2 v 1 {x i )) n . (80) 



We see that the integrals I n are the functionals only of Euler's variables. The 
consequences of this property will be discussed in the next section. 

In the case of 3-dimensional space we can construct the analogous integrals 
of motion, integrating the products of the vector (|68l) : 



K n , n .., n = J D <?tR il {ti)R h {ti):.Rj n (Si) (81) 

Changing variables ^ — > above we shall obtain: 

m) - mw) = ^Mn^))Am^) dp f x {x)) - (82) 



in 



In the 3-dimensional case the matrices AJ'{^{x)) satisfy the equation: 

e jkl A^(x))A^(x)) = e mnr ^^-detA(ax)), (83) 
Therefore (l82l takes the form: 

*,<{<*)) = ^ t „^^ (84) 
p(xj) ax r ox n 



"in 
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The integrals K n in (|8"T1) become 

K juj2 ... jn = Jj i xp(x l )R 3 Mx))RrMx))---R 3 Mx)) (85) 

Apparently, due to the presence of 9 g^ in (|84l) . we can not in this case 
express (l85l) only in terms of velocity and density, i.e. the Eulerian description 
does not admit this kind of integrals of motion. 

In the 3-dimensional case there is one more integral, which does not 
exist in any other dimension. Recall that the vector </&(£;), given by ( |67l ) 
is conserved modulo gradient, therefore the integral 

Q = f d 3 £J k &)R k &) (86) 

is conserved because R(£,i) = curl.J{^i). Transforming the ^-dependent variables 
to the £j-dependent ones in (l86l) we obtain helicity functional : 

Q = Ij^WvAm f^f' ' • (87) 

5 Euler variables. 

This section we shall discuss the relation of Lagrange and Euler description 
of hydrodynamics (see also a very deep and interesting discussion in [14J, 
there also could be found numerous references to the earlier investigations). 

In our approach, presented in Section 3 the phase space of the 3-dimensional 
fluid is 6-dimensional (xoo), as it naturally follows from Lagrange description. 
This could be compared with recent papers [IT], where (in our notations) only 
p(xi) and Vj(xi) are regarded as the coordinates in the phase space (this point 
of view of could be found also in different text books and articles) see for 
example [4], [9], p]. 

According to the conventional point of view the state of the fluid is 
determined by its velocity and density and therefore all other variables like 
ours x(C,i) are not needed. Indeed, solving the Euler equations of motion we 
can express the velocity v(xi,t) and the density p(xi,t) at time t via the 
initial data v(xi,0) and p(xi,0) . Then, from the definition of v(xi,t) Q 
follows 

*;(£(£,£),£)) = £(£,£). (88) 
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(here we have suppressed the indices of the arguments for brevity) . Apparently 
we can solve these equations with respect to x(£,t), provided we know the 
initial data a?(£,0). Therefore it may seem that the variables x(£, t) are 
unnecessary, as it could be obtained through the others. But the initial data 
x(£, 0) are half of the canonical variables in Hamiltonian formalism. So, in our 
approach we indeed need more variables. These additional variables provide 
the complete description of the fluid in a sense that solving the equations of 
motion we define not only the v(xi,t) and p(xi,t) but also the trajectories 
of the particles which could not be obtained in the conventional formalism. 
The situation is analogous to the rigid body rotation: here the phase space 
T is formed by 3 angles and 3 components of the angular momentum Jj. 
As the Hamiltonian depends on the components of Jj only, we can consider 
separately the evolution of the angular momentum. This is incomplete description 
for this mechanical system. The complete one certainly should includes the 
evolution of the angles, which define the location of the rigid body in space. 

In the case of fluid dynamics the complete description is to be done 
in the 6-dimensional phase space T formed by x(£) and p(£) (or £(x) and 
l(x)). If, as it is usually the case, the Hamiltonian depends only on l(x) and 
p(x), partial description in terms of the velocities and densities considered 
as "relevant "variables is possible. This means that we do not care about 
the evolution of the whole set of coordinates of the phase space which are 
considered as inessential. The "relevant "part of the coordinates does not 
necessarily form a simplectic subspace in T. The non-degenerate Poisson 
brackets in T could become degenerate on the subset of T, corresponding to 
the "relevant "variables. This is indeed the case in the rigid body and in the 
conventional fluid dynamics. In the first case the degeneracy of the algebra 
of Poisson brackets for the "relevant "variables — angular momentum is well 
known. Its center element is the Casimir operator of the rotation group. 

In the case of fluid dynamics the algebra of the "relevant "variables — 
velocities and densities for arbitrary dimension has the following form: 

1 

{vj(xi),v k {yi)} = ~ ^ (VjUfc(xi) - V k Vj(xi))5(x - y) 
{vj{xi),p{yi)} = —Vj5(x-y) 

{p(x t ),p( yi )} =0 (89) 

This algebra could be obtained from equations ( l3~Tl) . (|36l)and ( l37l ). The center 
of this algebra is infinite-dimensional and its structure depends on the dimension 
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of x-space. To find the center of it the following consideration could be 
useful. In the previous section we have considered the invariance of the fluid 
dynamics with respect to the infinite-dimensional group of special diffiomorphisms 
SDiff. This group acts on the Lagrange variables, and in [T4] it was called 
"relabeling symmetry". Apparently, this transformations do not affect the 
Euler's variables v(xi) and p(xi). Therefore these variables should have vanishing 
Poisson brackets with the integrals of motion, corresponding to the invariance 
with respect to the special diffiomorphisms. Indeed, the direct calculation 
gives: 

{R k i{£i),v m {xi)} =0, 

{R M (&,p(xi)} =0, (90) 
where we denote as Rki(£i) tensor 

Rki&) - (91) 

which contrary to the vector Rk(£i) defined in the previous section, is appropriate 
for any dimensions. Needless to say that the circulations (1731) also commute 
with v(xi) and p(xi). This observation gives us a guideline for the construction 
of the cental elements of the algebra ( 1591 . Indeed, if we could find a functional 
of Rki{£i), such that after transformation to x-dependent representation it 
will depend only of v(xi) and p(xi), then it automatically falls into the center 
of (1591 . As we have seen in the previous section, in case of 2-dimensional fluid 
we were lucky and the integrals I n , given by ( 1501 were expressed only via 
Euler variables. Therefore I n form the center, the fact which is known since 
long , see [4], [TJ] for the discussion. Here we want to note the remarkable 
fact: the symmetry, which is unknown to the variables v(xi) and p(xi), helps 
to find the center of its algebra! 

The 3-dimensional case we shall consider in details in the next Section. 
The Casimirs which follow from the reasonings, given above are the "helicity"functional 

Q = J Sxe jklVj (x)W k vi(x). (92) 

The other Casimir is the total number of particles iV(valid for any dimension, 
for d = 2,N = J ): 

N = f d 3 xp{xi) (93) 
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6 C 2 Hydrodynamics 



The 3-dimensional case, which is very important for applications was considered 
by many authors starting from XIX century. It is hardly possible to give an 
exhaustive list of references. Recently it was discussed in p2] (see also [H] 
for earlier references) where it was suggested to build Hamiltonian formalism 
for 3-dimensional Euler fluid using Clebsh parametrization for the velocity : 

v{ Xi , t) = da{ Xi , t) + P&y(xi, t) (94) 

where the new functions a(xi, t),/3(xi, t), 7(2^, t) together with density p(xi, t) 
are used for the construction of the coordinates of the phase space. What 
we are going to suggest here is an alternative approach, which has certain 
advantages [15] . 

Let us consider a mechanical system which is described by a pair of 
complex coordinates which belong to C 2 x 00: u a (xi,t), u a (xi,t), where 
a = 1,2. The Lagrangian for this system we shall take in the following form: 

L = J Sx{ l ^{u{xi,t)u(xi,t) - u{xi,t)u{xi,t)) 

(u(xi,t)du(xi,t) - du(xi,t)u(xi,t)) 2 + \\\ fn a\ 

+ m ^7 7^—, — 7^ K(w(x i ,t)M(x i ,t))},(95) 

QU{Xi, t)U{Xi, t) 

where we assume the summation over indexes. The canonical momenta, 
corresponding to the variables u a (xi), u a {xi) are given by equations (from 
now on we again will suppress the time argument) 

ui \ im - I \ 

Pa(xi) = -y««W- (96) 

As it expected for the Lagrangian which is a linear function of velocities, the 
equations (l96l) define the constraints on the canonical variables: 

A* (a*) = pl(xi) - —u(xi) a ~ 0, 

- 1 771 

= Pafa) + ~2^i)a ~ 0. (97) 

The Poisson brackets of the constraints are non-degenerate 

{A^(xi), XpiVi)} = im5 a p5(x - y). (98) 
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and we can use these constraints to eliminate canonical momenta p^(xj),p^(xj) 
using Dirac procedure [IB]. The resulting Poisson (Dirac) brackets for the rest 
of coordinates of the phase space T are: 

i 

{u a (xi),up(yi)} = —8 a p8(x - y), 

m 

{u a (xi),up(yi}} = 0, 

{u a (xi),up(yi)} = 0. (99) 

The Hamiltonian, corresponding to the Lagrangian ( 1951) has the following 
form 

f r (u(xi)du(xi) — du(xi)u(xi)) 2 T ,. . s . ... 

J ' QU[Xi)U[Xi) 

Now we shall explain why we consider this system. Let us form the 
following objects: 

v(xi) = —{u{xi)du{xi) - du{xi)u{xi)), 

p(xi) = u(xi)u(xi). (101) 

The notations we have used here are not accidental. The point is that if we 
shall calculate the Poisson brackets for (jlOip , using (j99|) the result will exactly 
coincide with l j89]h The Hamiltonian if, given by (11001) . being expressed via 
v(xi) and p(xi) will take the following form: 

H = J (fx^mpix^vixi) 2 + V(p(xi))}, (102) 

which also coincides with Hamiltonian given by (f47j) . 

The equations of motion for the variable i) have the usual 

form: 

Uct{Xi) = {H,U a (Xi)} 

u a {xi) = {H,u a (xi)} (103) 

Apparently, the correct equations of motion , including the continuity equation 
for variables v(xi) and p(xi) follow from ( 1103ft . 
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As was mentioned above, the description of the fluid in terms of u a (xi), u c 
is rather similar to the description which uses Clebsh parametrization. Indeed, 
these variables could be presented in the following form: 

Uafa) = sjpix^^ ( ^ {xi)/2s . n ^ ] (104) 



from where we obtain the representation for the velocity through angles 

4>(xi),ifj(xi) and a(xi) 

= 7^(d(f)(xi) - dip(xi)cosa(xi)) (105) 

These equation defines the velocity, if Clebsh parameters are known. Also, 
as is well-known (see e.g.[3j,[l2]) any differentiable vector field v(xi) has the 
local representation (11051) . In other words, knowing v(xi), we can construct 
Clebsh parameters a(xi),<p(xi),ip(xi) with some ambiguity. This ambiguity 
arises as a set of integration constants . In our construction this ambiguity 
could be understand as follows. The Lagrangian function ( |95i l we consider is 
invariant with respect to the symmetry group U(2) which acts as follows: 

u a (xi) -> u a (xi) = T al3 up(xi), T + T = 1 (106) 

and according to Noether's theorem the integrals of motion, which is the 
generators of these transformations are : 

t° = J d 3 x^u(xi)u(xi), t a = J d 3 xu(xi) Y u ( x i) ( 107 ) 

The transformations 01070 change the Clebsh variables, but does not affect 
the Euler's variables. So, in particular, the constant shift of the angle <f>(xi) 
is generated by ex-Casimir N, which in T has lost its status, the generator 
t 3 shifts the angle ip(xi) , the other two generators mix the angles ij)(xi) and 
ot{xi) . So the system described by variables (u a (xi), u a (xi)) is a Hamiltonian 
system with symmetry and we can reduce its phase space by procedure given 
by Souriau p6j and Marsden and Weinstain [T7] . The reduced phase space is 
the space, where "live"almost all the Euler variables. The latter means that 
the procedure of reduction implies fixing the integrals of motion, in particular 
t° = y does not anymore belongs to the set of variables. 

The only problem we have now is the "helicity"functional, which still the 
is the Casimir and the reduction does not remove it. For finite dimensional 
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systems the existence of Casimir implies the degerancy of Poisson brackets. 
It could be easily seen from the following consideration. By definition the 
Casimir C should has a vanishing brackets with all variable 

{p k ,C} = 0, {q k ,C} = (108) 

where pk,qk are all set of coordinates of the phase space. If the Poisson 
brackets are non-degenerate, the equations (11081) mean that C is a constant. 
The situation for the infinite dimensional system is different because of 
existence of so called the functionals with zero variation . Consider for example 
an infinite dimensional system , which is described by the set of canonical 
variables p(x),q(x) where x G R . The Poisson brackets are non-degenerate: 

{p(x),q(y)} = 5(x-y). (109) 

In this case we can easily construct a nontrivial functional, which will have 
vanishing Poisson brackets with all variables p(x),q(x). It has the following 
form: 

r oo p'{x)q{x) — p(x)q'(x) 



C = dx ' ^gpi-i^g ^ (no) 
p 2 (x) + q 2 (x) 



and has a meaning of a winding number for the phase of the complex variable 
a(x) = p(x) + iq(x), i.e. C is what physicists used to call topological charge. 
Note that in order for C to be the Casimir it is not necessary to impose the 
condition on a(x)| I ^_ 00 = a(x)\ x -> +oc and compactify R. In this case C will 
take an integer values and indeed will be the winding number. 

The "helicity"functional has the same origin as the functional C in this 
example. In order to see it let us introduce a unit four vector F k 



1 ( Uxixi) \_( Fi(xi) + iF 2 (xi 



(111) 



This four vector maps S 3 — > S 3 , provided we impose on the variables u a (xi), u a { 
the asymptotic conditions : u a (x) — > u° a , when \x\ — > oo and compactify R? . 
The helicity functional Q given by ( 1921) could be written in the following 
form: 

Q = g J ^xeabcdeijkFadiFbdjFcdkFd, (112) 

which is the standard representation for the winding number of the map 
S 3 — ► S 3 , so called Hopf invariant. Here again we should note that even if we 
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neglect the asymptotic conditions on u a (x) together with compactification of 
.R 3 , Q is invariant with respect to local variations and therefore has vanishing 
Poisson brackets with u a (xi), u a (xi). So, the conclusion of these arguments 
is that for infinite dimensional mechanical systems the existence of Casimirs 
does not necessary implies the degeneracy of Poisson brackets, provided these 
Casimirs are related to the geometric properties of the phase space and the 
helicity belongs to this class of "friendly"Casimirs. 



7 Inclusion of the electromagnetic interaction. 
Plasma. 

We shall consider plasma as a fluid of two components, namely electrons with 
mass m and electric charge (— e) and ions with mass M and charge (+e). 
The coordinates of electrons we shall denote as x(^), while the coordinates 
of ions will be The interaction of the components of the plasma with 

the electromagnetic field A^x) is governed by the following Lagrangian: 



L — Lg + L i Q n + ejd 3 £ 



-ejd 3 E 



-\ J d 3 xF^(x,t)F^(x,t), 
where Lq 1 and L™ n are "free "Lagrangian 



(113) 



Lt = J d 3 C 
}™ = J d 3 E 



mx fat) 



MX(E h t) 



f e \det 



- f lon (det 



dxj(£j,t). 
dXj(Ei,t) 



2 ' dE k 

and F^ u (x,t) denotes the electromagnetic field tensor: 
F^(x, t) = d^A u (x, t) - d^A^x, t) 



(114) 



(115) 



The advantage of Lagrangian description of fluid (plasma) is clear. Using of 
the coordinates of charged particles as the fundamental variables allows to 
introduce the interaction with electromagnetic field. 
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The Lagrangian (11131) possesses U(l) gauge invariance. We shall use the 
usual Hamiltonian formalism for the constraint system [18J. The canonical 
variables are: 

6L 





m) = - 
< 


*(S,), 


P&) = 


A{ Xi ), 


Pem {.Xi) 


A {xi), 


Pem{ X i) 



SL 



SX(l 



5L 



SA(xi) 



-E(xi 



0. 



(116) 

where E{x r ) is the electric field strength: 

E{x l ) = -VA {x i )-A{x i ). (117) 

The last of the equations (I116P is actually the primary constraint. The 
Legendre transformation of the Lagrangian ( 11131 ) gives us the canonical 
Hamiltonian: 



H = fd 3 x[± (Pl n (xi) + H\xi)) + A { Xi )VP e , 



I d ^[t (Mi) + e^(x(6))) 2 - eA (x&)) + f^det^ 



eA(X(EA)) + eA {X{E l )) + f™(det a -^g^] 



(H8) 



where H{xi) = curlA(xi) is the magnetic field strength. The requirement of 
the conservation of the primary constraint P e ° m (xj) = gives the secondary 
constraint ( Gauss law ): 



Vj {P e m{Xi))j ~ e (pel{Xi) ~ Pion{xi)) 

Now we can add to the primary constraint P® n 



0. 



(119) 



the gauge fixing 



condition A {xi) = and eliminate these variables from consideration. Introducing 
the x-dependent functions instead of £ and S-dependent, as we did in the 
2-nd section we obtain the Hamiltonian of the plasma in the following form: 

(l(xi) + ep e i(xi)A(xi 



H = Jd 3 x[l (P*Jxi) + H^Xi)) + 
(L(xi) - ep ion (xi)A(xi) 



2mp e i(xi) 



2mpi 



+ V e l(p e l(Xi)) + V ion (p ion (Xi)) 



(120) 
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where we have introduced 



L{ Xi ) = fd 3 ~P(~ l )6(x-X(~ i )) (121) 

This Hamiltonian is gauge invariant with respect to the transformations, 
generated by the constraint (|119p . Further we impose the Coulomb gauge 
condition on the electromagnetic field 

V^A,^) = (122) 

and following the usual procedure will eliminate the longitudinal components 
of A(xi) and P em (xi). As a result of the gauge fixing, the longitudinal part 
of P em (xi) give rise to the Coulomb term in the Hamiltonian 

H C ol = -e 2 / d 3 X {p e l{Xi) - p ion (Xi)) (-——; ~A {peliVi) ~ PioniVi)) ■ 

J \ Att\x — y\ J 

(123) 

and the whole Hamiltonian takes the following form: 

H = Hcoi + I ($x[\ (P? m± (xi) + H\x 

2mp e i{xi) 

(L(xi) - ep ion (xi)A ± (xi)) , 

+ ^ = r~\ + V ion(Pion(Xi)) , 124 

where the subscript J_ denotes the transverse components of electromagnetic 
variables, for which the Poisson (Dirac) brackets are given by [18] 

{PlJ(x), A ±k (y)} = (s jk - 5(x-y). (125) 

In equilibrium plasma the Coulomb interaction is known to be screened by 
the cloud and the residual Debye interaction is a short rang one (see e.g. [i~9]). 
In [20] it was stated that Coulomb term is reduced to the local functional of 
the difference of charge densities (p e i(xi) — Pi on ( x i))- 

Following our consideration of a fluid in the Section 3, we can introduce 
the canonical coordinates £(xi),jr(xi) for electron and E(xi),U(xi) for ion 
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components of the plasma: 



{nj(xi),€k(yi)} = $jk5(x-y), 
{Uj(xi), ^kiVi)} = djk$(x - y) 



(126) 



for which 




fa) 



d7r(xi) 




dxj 



(127) 



dxj 



Substituting (11271) into (11241) we obtain the Hamiltonian in terms of the 
canonical variables. 

The Lagrangian (11131) is invariant with respect to the volume preserving 
diffeomorphisms of both components of plasma separately, so we shall have in 
this case two sets of conservation laws - one for electrons, the other for ions. 
Applying the procedure, which we have described in the previous section we 
shall construct the conserved circulations: 



where lj(x) and Lj(x) are given by (1127ft . In the same way we can construct 
the analogues of the integrals (l85l) and (1871) for this case. Note that in the 
case of plasma the equation (1441) is not valid due to the presence of the 
electromagnetic field, instead we have: 



therefore the equations (11281) are not the circulations of the velocities. In the 
same time, adding circulations V£ l and V™ 71 on the common contour A we 
obtain: 



so the total circulation of elections and ions velocities is conserved even in 
the presence of the electromagnetic interaction. 




(128) 




(129) 




(130) 
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Summarizing we can say that the phase space of plasma T in Coulomb 
gauge is the space with coordinates A±(x), F^ m (x); 7?(xj); Tl(xi) 
with Poisson brackets given by (11251) and (I126P . The evolution of a state in T 
is defined by the Hamiltonian (I124p . In the same time we can use incomplete 
description in the terms of Euler variables, but in the case of plasma the 
Poisson brackets for the velocities are changed because of the presence of 
electromagnetic field in equations (11291) . For example the Poisson brackets 
for the velocities of electrons take the following form: 

djV k {xi) - d k Vj{xi) + MdjA k {x) - d k Aj{x)) 5{x-y). (131) 



mp(x) 

The same way are changed the Poisson brackets of the velocities of ions, while 
the brackets involving densities are not changed. 



8 Gravitating gas 

In the recent paper [6| we have considered the Hamiltonian formalism for 
fluid and gas based on the Lagrangian description. It was pointed out in this 
paper that apart from other advantages, the Lagrangian description, which 
uses the trajectories of the particles of fluid (gas) as the dynamical variables, 
is the most convenient for the introduction of interaction. In particular in the 
previous section it was demonstrated how the introduction of the electromagnetic 
interaction of particles which constitute the fluid, provide us with the theory 
of plasma. In the present paper we are going to consider in analogous way the 
theory of gas of particles which interact with each other through gravitational 
Newton potential. The system of such particles could be considered as a 
model for the motion of stars in a galaxy when the gravitation interaction 
prevails all other interaction. The total number of stars in typical galaxy is 
of the order of 10 13 — 10 14 , so it may be reasonable to consider this collection 
of "particles"as a gas. 

The simplest model, which is usually used for numerical simulation of 
N-body model of galaxy is described by the Hamiltonian 

N -.2 N 1 

i=l 2 m i i^j \ x i x j\ 

where Xi,pi are canonical coordinates of particles (stars) with masses rrii. 
The model we are going to consider is based on the dynamics described by 
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H with the assumption = m, when iV — > oo. For this limit there appears 
a natural desire to consider a continuous distribution of the particles as it is 
done in the theory of fluid or gas [22j. 

The Lagrangian of the continuous system of particles interacting through 
Kepler potential has the following form: 



where 7 denotes the gravitational constant. Comparing |132| with [19] we 
see that the only difference is in the form of "potential energy"therefore the 
canonical formalism for gravitating gas essentially the same , as for fluid. In 
such a way we can immediately write the Hamiltonian of gravitating gas: 

J ^ 2m 2 J s s |x(6,t) -x{&t)\ V ' 

Inserting into the integrals in the r.h.s of |133| the unity 

1 = J d 3 xS(x- x^i)), (134) 

as we did in the Section 3 we can rewrite |133| via x-dependent variables: 

Hr-\ 7 f j3_j3..P(xi)p(xi 



H = / d s x- —A^Xi) - ± / d 6 xd s y ^\:^\ lJ . (135) 

J 2mp{Xi) 2 J \x — y\ 

Certainly the attractive interaction between particles will dramatically change 
the properties of gravitating gas 

The equations of motion, which follow from the Lagrangian (11321) have 
the form: 

mSfc, *)+7/ *fMS) |ggj;g, - (136) 

Translating equation ( 11361) on to the language of Euler variables, as has been 
done above we arrive at the following set, including the continuity equation: 

ff(x i} t) + v k (xi,t)^-^x i} t) = ^^Id'y^^l (137) 
oxk max J I x — y \ 

d 

p(xi,t) + (p(xi,t)v(xi,tf) = 0. (138) 
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Note that the r.h.s of the equation (16) in the case of ordinary gas or fluid 
is expressed through the internal pressure p(xi); 



V{Xi, t) + V k [Xi, t)- V{Xi, t) = — ——p{Xi) 

OX k p(Xi) ox 



(139) 



The set (16), (17) defines the evolution of initial distribution of p(xi, t ), v(xi, t ) 
of gravitating gas and, besides it could be used to find the static configuration 
of this gas for different boundary conditions. In particular we can explore the 
possibility of the existence of the static isolated configurations of gravitating 
gas. Isolation here means that density p(xi) vanishes at infinity. Note, that for 
usual fluid, gas or plasma these kind of solutions are forbidden due to virial 
arguments, known in the case of plasma as Shafranov's theorem [23].|24|. For 
the gas, describing by equation (11391) this theorem could be proven as follows. 
First, using the continuity equation let us rewrite equation ( 11391 ) for static 
case in the following form: 



_d_ 

dxi 



p(xi)v k (x i )v j (x j ) + S jk p(xi 



0. 



Integrating (11401) with Xj over R 3 we obtain: 





= Jd 3 



xx 



3 dx k 



Jd 3 







x 



dx k 
J d 3 x5jk 



p{x i )v k {xi)v j {x j ) + S jk p(x 
p(xi)v k (x i )v j (x j ) + 5 jk p(xi 



(140) 



(141) 



p(x i )v k (x i )v j (x j ) + S jk p(x 



For usual gases p{xi) ~ p 7 (xj),7 > 0, therefore the integral over divergence 
will vanish for isolated solutions for which p(x{) — > when \x\ — > oo and we 
the obtain the following equation: 



d 3 x 



p(xi)v 2 (xj) + 3p(xi) 



0. 



(142) 



Apparently, this equation could be satisfied only for the case p{xi) = 0, i.e. 
there is no isolated in the above formulated sense, static solutions of the 
equation (11391) . In order to have a static solution of (11391) we need to change 
the boundary condition p(xi) — > to the condition p(xi) p as , where p as 
is asymptotic uniform density. 
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Now we shall show that the arguments of this theorem bring no obstacles 
for gravitating gas. For this we again will rewrite the equations (16) for the 
static case, using continuity equation in the following form: 



d 

d Xj 



p{Xi)v k (Xi)Vj(Xi) 



1 ( X 9 

-p{Xi 



m dxk 

Integrating (11431) with x k over R 3 we obtain: 



d A y 



\x-y\ 



0, 



(143) 



d oc 



_d_ 

dxj 



1 

p(x i )v k (x i )v j (x i ) p(x 



d 



m dxi 



d 6 y 



\x - y\ 



0. (144) 



Consider the first term of the integrand in (11441) . Integrating by parts as 
above and taking into account the asymptotic conditions for p(xi) we obtain: 

d 3 xx k -^(p(xi)v k (x i )v j (xi) ) J = - J d 3 xp(xi)v 2 (xi). (145) 

Integration of the second term of the integrand in ( 11441 ) is straightforward, 
yielding 



m 



_d_ 

dxi 



3„.p(yi,t) _ 7 



d 6 y 



x — y\ 2m 



73„J3„.P(^)P(^ 



d xd y- 



\x-y\ 



(146) 



In such a way from equation fll44[) we obtain: 



d 3 xp(xi)v 2 (xi) — I d 3 xd 3 y 



p(Xi)p(Xj 

\x-y\ 



0. 



(147) 



This relation apparently could be satisfied for non-trivial configurations of 
p(xi),v(xi). The energy functional, corresponding to the Lagrangian ( IPTl ) has 
the following form: 



in 



1 



E = — d 3 xp(xi,t)v 2 (xi,t) - d 6 xd 6 y 



p(X j,t)p(yi,t) 
\x-y\ 



(148) 



The two terms of the energy functional have clear interpretation as kinetic T 
and potential U parts of energy and equation (11441) expresses famous "virial 
theorem" p]: 

2T= -U. (149) 
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Note here, that in the "virial theorem "equation (1149p holds true for mean 
values of kinetic and potential energies, while in our case of static solutions 
there is no need to average over time. 

Using the relation (11491) we can easily find the total energy of static 
configuration of the gravitating gas: 

E static = ~— I d 3 xp(xi)v 2 (xi) = (150) 



2 

7 f jh p(xi,t)p{yi,t) 



J d 3 xd 3 y 



4 i ' \x — y\ 

So, the total energy of the static solution is negative, as for the bound state 
of Kepler problem. 



9 Properties of Static Solutions 

The equations which define our static configurations of gravitating gas have 
the following form: 

Vk{Xi)^ — v{Xi) = — — / d y- 



ctefc max J \x — y\ 

-^L(p(x t )v(x t )) = (151) 

Taking divergence of the first equation : 



d(v k (xi)—v(xij) — —A / d 3 y 



dxk ' m J \x — y\ 

= l(-4n)p( Xi ), (152) 
m 

we can write the whole set of the equations for static configuration in a pure 
local form: 

d[v k (xi)- — v(xi)} = -4ir—p(xi), 
K oxk ' rn 

d 

dx (y k (xi)-^-v(xi)) = 0, (153) 
31 



where the last equation requires the expression v k {xi)-^-v{xi) to be a gradient. 

Now we are going to derive an important inequality, which bounds the 
{—E static) °f an Y solution of (1153P from below. For this let us introduce the 
notation for the potential U(xi): 

U( Xi ) = - I d 3 y-$^- (154) 
J \x — y\ 

Integrating by parts we obtain the following relation: 

/ d 3 x(dU( Xl )) 2 = 4vr / d 3 xd V {X \ , ^ i,t) (155) 

Using (I155H we can write the expression for the energy given by (11501 ) of any 
static configuration, which is the solution of (11531) in the following form: 

= -j^ J d 3 x(dU(x t )f (156) 

Substituting into (11560 the expression for the gradient of the potential U(xi) 
from the first equation (11511) . we obtain the expression for the E static only 
through the field of velocity: 

- E static = J d 3 x(v k (xi)-^-v(xi)y (157) 

This form is most convenient for the derivation of desired inequality. Now let 
us consider a function f(xi) from the Hilbert space which consists of all 
measurable functions on R 3 , which have at least one derivative and square 
integrable on R 3 together with its derivatives. In particular we assume that 
the density p(x^) belongs to W\. Taking into account the equation (11521) we 
have 

\ J d 3 xp(x i )f(x l )\ = \ f d 3 xf(x i )d k (v j (xi)d j v k ( Xi 
= li^ / d 3 xd k f{x i )v j {x i )d j v k {x i )\ 

^(jd'xidkfix^^Jd'xivjix^dM^)) 2 ) 112 , (158) 



< 



where on the last step we have used Cauchy inequality. From f|158f) we 
immediately obtain the inequality: 

(/ d 3 xp(xi)f(xi)) 2 

~ aHc ^ *7 fd 3 x(d t f(Xi)) 2 • <159) 
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which is valid for any f(xi) from W\ and is saturated for f{xj) = U(xi). 
Indeed, let us calculate the derivative of the functional in the r.h.s of (11591) 
with respect to f(xi): 

5 (l <Pxp(xi)f(xi 



Sf(xi) Sd?x{d k f{ Xi )y 

, 2 



2 



(j<Pxp(xi)f(xi)) r p (x.) Af( Xi 



+ 



(160) 



Sd?x{d k f{xi)Y \.(I&xp{xi)f{xi) Jd3x(d k f(x t ))i- 
This derivative vanishes for f(xi) = f ma x( x i) given by 

fmaxiXt) = CA-V(Xi) = C'U{Xi), (161) 

where C and C are inessential constants, which do not enter into the functional. 
It is easy to prove that the second variation of this functional is negative 
on the fmax(xi), so this function provides the absolute maximum for the 
functional and due to the equations (1155H157I) its value coincides with the 
l.h.s. of ( 11591 ). However, the function f ma x(xi) does not belong to the Hilbert 
space W2, because the potential U(xi) has the asymptotic behavior |4 at 
infinity and therefore is not square integrable. 

The integral J which enters into inequality ( 11591) could be written in the 
following form 

J= f ^xpix^Uixi), (162) 



where we denoted as /jv(a?i) the normalized function f(xi): 

yJfd?x(d k f(Xi)y \\dkf(Xi)\\2 

In the Hilbert space W\ the Holder's inequality holds true: 

d z xf{x i )g{x l )\ < \\f\\ v \\g\\p>i - + ^ = 1 (164) 

as well, as the remarkable Ladyjenskaya's [25] inequality: 

< (48) 1 / 6 ||a fe /|| 2 (165) 
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Here we use the standard notations: 



11/11, = (jd*xf{ Xi y) llP , 

\\dkf\\ P = (J d 3 x\d k f(x^) 



i/p 



From these two inequalities we obtain the following bound for the integral J: 



The facts mentioned above could possibly support the statement that the 
(—E static ) is bounded from below by appropriate norm of the density. Indeed, 
the absolute maximum of the functional | J\ should be bigger when its maximum 
in a restricted space like given by (11661) . However we can not present the 
rigorous proof of this statement. 

One of the other general property of a static configurations of gas or fluid 
is the existence of the topological charge — "helicity"(or Hopf invariant), 
which explicit form is 



This object is not only the integral of motion of the equations (16) but it also 
is the central element of the algebra of Poisson brackets of [v[x,i), p(xi)) [6]. 
The existence of such an object brings additional argument for the stability 
of the solitons. The role of "helicity"in the case of the solitons in plasma 
was pointed out in [2T],[24]. Moreover, in [?] it was shown that there exists 
a remarkable inequality which bounds the energy of plasma solitons from 
below by q 3 ^. The derivation of such inequality for our case (and in general 
for fluid solitons ) is highly desirable and we going to consider this question 
in the future publications. 

10 On the Possible Structure of the Static Solutions 

The equations (I153f) which define the static configurations of gravitating gas 
are 3-dimensional nonlinear partial differential equations and the probability 
to find an analytic solution is very low. The only case where a class of 
solutions was found in a similar situation is the t'Hooft-Polaykov monopole, 
but there the requirement of the spherical symmetry simplified essentially the 
problem. In our case we can not expect the spherically symmetric solution 



A < (48) 1 / 6 ||p|| 6/5 



(166) 




(167) 
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because the continuity equation requires the trajectories of the particles 
be closed in order to provide the static configuration for v(xi),p(xi). The 
simplest and most symmetric configuration we could expect for our case is 
the toroidal structure, where the density is concentrated in the vicinity of the 
axis of the toroid, while the field of velocity is tangential to the embedded 
one into the other toroidal surfaces. It is not the first time the toroidal-shape 
soliton appeared in the context of the theory of continuous media. Since 
the pioneer works of Lord Kelvin in XIX century to the present time it was 
studied by many scientists both mathematician and physicists and recently 
the interest to the subject was again attracted by the works of Faddeev and 
Niemi [20], [21] 

It is clear that in the case of attractive gravitational interaction the 
particles, which constitute the gas move around the region with bigger density 
(like planets around the Sun) so there should exists a collective motion in the 
approximation when the radius of torus tends to infinity and we can speak 
about cylindrical rather when toroidal configuration. Indeed, let us consider 
this axially symmetric tornado-like solution of ( I153p . For that we shall write 
the Ansazt: 

v(xi) = v(r)(--,-,0), 

p( Xi ) = p(r), (168) 

where r = ^/x 1 + y' 2 . The second and the third equations (11531 ) are satisfied 
by (11681) . while the first gives 

1 T 

-d r v(r) 2 = 4tt— p(r). (169) 
r m 

This solution shows that the density stays an arbitrary function (what is 
expected for the partial differential equations) and the velocity grows with 
radius up to its asymptotic value. The later is the consequence of the approximation — 
here we actually have 2-dimensional potential logr instead of Coulomb K This 
example is to demonstrate that the particles which constitute the gravitating 
gas in their collective motion form the localized object. 

The traditional way to tackle 3-dimensional gas or fluid is to introduce for 
the velocity Clebsh parametrization suggested in [8] and recently discussed 
in|ll|. In general case this parametrization has the following form: 

v( Xi ) = f{xi)Sg{xi) + dh(xi), (170) 
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where f(xi), g{xi), h{xj) are scalar functions. For the toroidal solution with 
z - as the axis of symmetry we shall assume that the density p{xi) does 
not depend upon the azimuth angle <f> and Clebsh parametrization takes the 
form: 

v(xi) = cosa(r, z)df3{r ) z) + kd<p(xi), (171) 
where r, z, (p are cylindrical coordinates, 

r = \J x 2 + y 2 

4>(xi) = arctang— 
x 

In fluid dynamics the parametrization (1171H has an interesting mechanical 
interpretation which we shall discuss elsewhere. The "helicity"functional for 
this case has the following form 

q = J kd<p(xi) A dcosa(r, z) A d/3(r, z) 
= 2nk j dcosa(r, z) A d(3(r, z) (172) 

The function f3(r, z) in (11711) should has a singularity at the point (r c , 0) in 
the semi-plane (r, z) where r c is the radius of axis line of the torus. When 
the point (r, z) goes around (r c , 0), the function (3(r, z) increases on 2ti0. The 
gradient of such function will be the vector tangential to the toroidal surfaces. 
The addition of the gradient of azimuthal angle (f> makes the velocity (11711) 
winding also in the azimuth direction. Such a behavior of the field of velocity 
leads to the nontrivial Hopf invariant, which characterizes the homotopy 
classes 7Ts(S 2 ). Described this way the Clebsh parameters, together with the 
density p(xi) could be found by numerical integration of the equations (11531) . 

The construction we presented and discussed above would be incomplete 
and too academic without an example of the galaxy which indeed may 
demonstrate the toroidal structure. Fortunately such galaxy does exists. It 
was discovered back in 1950 by Art Hoag and recently a very good picture 
was obtained by Hubble telescope (see figure 1). On this picture clearly seen 
the toroidal structure. In the gallery of the galaxies which is available on the 
sites |http:/ /ww w.astro nomy.com| or |http: / /hubblesite.ogr we can find some 
other examples with the form more or less close to the torus, but the Hoag's 
object is the best of all. 

1 Compare this with the similar parametrization of the dynamical variables in [21] for 
the case of plasma. 
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Phc. 1: The above photo taken by the Hubble Space Telescope in July 
2001 reveals unprecedented details of Hoag's Object and may yield a better 
understanding. Hoag's Object spans about 120,000 light years and lies about 
600 million light years away toward the constellation of Serpens. 
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11 Concluding remarks 



Here we have considered different aspects of Hamiltonian formalism for several 
types of continuous media - gas, fluid, plasma and gravitational gas. In 
spite of different properties of these media, all of it could be treated within 
universal formalism and difference arises due to specific interaction between 
constituents. As it was shown the Lagrangian approach provides the simplest 
way to introduce different types if interaction because it uses coordinates of 
particles as basic variables. Having constructed Lagrangian formalism we can 
derive Euler equations of motion and also the whole set of Poisson brackets 
for Euler variables. 

The C 2 formalism of continuous media considered in Section 6 gives a 
frameworks for quantum theory which may be useful for the description of 
Bose-Einstein condensation in gas and this subject needs further in-depth 
study. 

Another subject which we are going to explore in the future is the theory 
of gravitating gas. Real challenge here is to find analytically toroidal soliton 
which could describe ring galaxies. 
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